HIGH RESOLUTION A=1MM CARMA OBSERVATIONS OF 
LARGE MOLECULES IN ORION-KL 
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ABSTRACT 

We present high resolution, Combined Array for Research in Millimeter- Wave 
Astronomy (CARMA), A=lmm observations of several molecular species toward 
Orion-KL. These are the highest spatial and spectral resolution 1mm obser- 
vations of these molecules to date. Our observations show that ethyl cyanide 
[C2H5CN] and vinyl cyanide [C2H3CN] originate from multiple cores near the 
Orion hot core and IRc7. Additionally we show that dimethyl ether [(CELj^O] 
and methyl formate [HCOOCH3] originate from IRc5 and IRc6 and that ace- 
tone [(CH 3 ) 2 CO] originates only from areas where both N-bearing and O-bearing 
species are present. 



Subject headings: astrochemistry- 
radio lines:ISM 



-ISM:individual(Orion-KL) — ISM:molecules- 



1. INTRODUCTION 

The Orion-KL region is the closest (~480 pc) site of massive star formation to us 
( jGenzel et al.lll98ll ). There are several cloud components (e.g. hot core, compact ridge, ex- 
tended ridge, and plateau) w hich are associated with Orion-KL, each with different chemical 
and physical properties (e.g. iBlake et al.lll987l ). The two most chemically interesting com- 
ponents, the hot core and compact ridge, are separated by only ~5800 AU. Large oxygen 
bearing species (e.g. methyl formate [HCOO CHc;] and dim ethyl e ther [(CH 3 ) 2 0]) have been 
observed primarily toward the compact ridge iLiu et al.l (e.g. l2002l ). while large ni trogen bear- 



ing s pecies (e.g. ethyl cyanide [C 2 H 5 CN]) are located toward the hot core (e.g. 
1987T) . 



Blake et al. 



The chemical pathways to most large detected interstellar molecules in hot molecular 
cores are not known or are poorly understood (e.g lOuan fc Herbstl 120071 ) . The spatial lo- 
cation of a molecular species, in relation to continuum, shocks, and other species, can give 
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indications of its formation mechanism. In order to further our understanding of these forma- 
tion mechanisms we have conducted high resolution observations of several large molecular 
species: acetone [(CH 3 ) 2 CO], dimethyl ether [(CH 3 ) 2 0], methyl formate [HCOOCH 3 ], ethyl 
cyanide [C 2 H 5 CN], methanol [CH 3 OH], formic acid [HCOOH], acetaldehyde [CH 3 CHO], and 
ethanol [C 2 H 5 OH]. 



Whi 



Liu et al. 



e many of these molecul ar species have been observed at moderate resolution (e.g. 



2002 



by iFriedel et al. 



Bla ke et al.lll996l ) (CH 3 ) 2 CO has not. Its original detection in this source 
(120051 ) was with a rather large beam (~ 12" x 9") and indicated that its 
emission was not associated with either the hot core or compact ridge. In this paper we 
present the highest spectral resolution 1mm observations of these molecular species to date. 



2. OBSERVATIONS 

The observations were taken in 2007 March with the Combined Array for Research in 
Millimeter- Wave Astronomy (CARMA). All data were acquired in two ~5 hour long tracks 
(one for (CH 3 ) 2 0, C 2 H 5 CN, (CH 3 ) 2 CO, and continuum and one for HCOOCH 3 , CH 3 OH, 
CH 3 CHO, C 2 H 5 OH, and HCOOH) with 9 antennas and a phase center of a(J2000) = 
05 ft 35 m 14 s .35 and 5(J2000) = -05°22'35".0. The typical synthesized beam was ~ 2.5" x 0.85 
(54 times higher resolution than the original acetone detection). The u — v coverage of the 
observations gives projected baselines of 5.8-196.5 kA (7.5-255.5 m). Thus any structure 
larger than ~ 16" will be r esolved out by the a rray. We compared our flux values to those of 



single dish observations of ISutton et al.l (119851 ) and determined that we are resolving out lit- 
tle flux (< 15%) from most of the lines we observed. The only exception is the J=8_i j g — 7qj 
transition of CH 3 OH, where we are resolving out ~90% of the flux. 

The correlator for the first track was configured to have two 500 MHz wide windows for 
continuum and four 32 MHz windows for spectral lines (half in each sideband), while the 
correlator for the second track was configured to have six 32 MHz windows (three in each 
sideband). The spectral line windows had 63 channels with a spacing of 488 kHz (~0.64 km 
s _1 ). Uranus was used as the flux density calibrator and 0530+135 was used to calibrate the 
antenna based gains. The absolute amplitude calibration of 0530+135 from the flux density 
calibrator is accurate to within ~20%. The internal noise source and 0530+135 were used 
to correct the passbands of each window. The da ta were calibrated , continuum subtracted, 
and imaged using the MIRIAD software package (jSault et al. 1995 ). 
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3. RESULTS 

Table [I] lists all the molecules that we searched for along with their parameters. The 
first column gives the molecular formula, the second column gives the transition quantum 
numbers, the third column gives the transition frequency in MHz, the fourth column gives the 
upper state energy (E u ) in K, the fifth column gives the line strength multiplied by the square 
of the relevant dipole moment {Sfi 2 ) in Debye 2 , the sixth column gives the synthesized beam 
size from the observations in arcseconds, the seventh column gives the 1 a rms noise level 
from the associated spectral window in mJy bm -1 , the eighth col umn given the perce ntage 



of missing flux when compared to the single dish observations of ISutton et al.l (119851 ). and 
the ninth column gives the literature reference for the molecular parameters. The table has 
been broken into 2 parts. The upper section lists the transitions detected in narrowband 
widows while the lower section lists transitions detected in wideband windows. 



3.1. Continuum 



We used two wideband (~500 MHz) windows to map the continuum of Orion-KL. 
Several molecular lines were detected in these windows (see §3.8p and were removed before 
producing the final continuum image. Figure [1] shows the 231 GHz continuum of Orion-KL. 
Known bright infrared and radio sourceaj are denoted with a "+" and are labeled. All maps 
in this paper are contoured in ±3 a, ±5 a, etc. contours. The 1 a rms noise level of the 
continuum is 11.3 mJy bm" 1 (1 a rms noise levels for the other maps are listed in Tabled]). 
Most of the continuum comes from the hot core (865 mJy bm -1 peak flux), however there is 
notable emission from IRc6 (232 mJy bm -1 peak flux), and there is also emission associated 
with IRc5 (254 mJy bm" 1 peak flux) and IRc7 (101 mJy bm" 1 peak flux). No emission 
was detected from either the compact ridge or BN, above our 1 a rms noise level. There is, 
however, a notable hole in the continuum around the location of BN. 



Previous 1mm array observations of this region by iBlake et al.l (119961 ) indicated that 
the hot c ore was compo s ed of at least 2 clumps and high resolution observations of ammonia 
[NH 3 ] bv lMigenes et al.l (1989) indi cate numerous clumps along the long axis of the hot core. 
Additionally, iBeuther et al.l (120041 ) detected several continuum peaks near the hot core at 
865 /im, which we were not able to resol ve. The elongation o f the continuum is also closely 
aligned with known SiO outflows seen by lWright et al.l (119951 ). Note that the emission from 



1 Some sources are not labeled as they are too close together, namely IRc2 and radio source I, near the 
hot core. 
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IRc6 is not a point source but rather a rough "C" shape around BN. There are also notable 
emission peaks both to the southwest and west of IRc5 which are not associated with known 
infrared sources. 



3.2. Methanol [CH 3 OH] 

Two transitions of CH 3 OH were detected in our narrowband windows. Figured shows 
the channel map of the t>LSR=5.5 km s -1 emission of the lower energy (J=8_i 5 s — 7qj) 
transition of CH3OH. The CH3OH map indicates some wide spread emission but has distinct 
concentrations near the hot core, IRc5, IRc7, and IRc6 (some are not seen in the map as they 
are at different velocities), but possibly most notable is lack of emission from the compact 
ridge. There is also a notable emission peak extending to the southwest from the hot core 
(hereafter hot core-SW) that corresponds to weak continu um emission. Convolving our 



results with a 29"beam and comparing them to those from ISutton et al.l (119851 ) indicates 



that we are resolving out a large percentage of flux ~90% from this transition. 

Figures [2b-e show the spectra toward IRc5, IRc6, the hot core, and IRc7 respectively. 
The highly asymmetric line profile of the IRc5 spectra indicates that CH3OH may be self 
absorbing near 9 km s -1 along this line of sight, while the IRc6 spectra peak at the same 
velocity as the IRc5 absorption. Both the hot core and IRc7 spectra show possible absorption 
features near 7.6 km s _1 . 

Fits to the spectral lines are given in Table [2j The first column lists the molecule and 
the second column lists the source. The third column gives the figure number for the spectra. 
The fourth, fifth, and sixth columns list the fit parameters: intensity in Jy bm _1 , Wlsr i n km 
s _1 , and FWHM in km s _1 , respectively. The seventh and ninth columns list the calculated 
opacity at the line center for assumed lower and upper rotation temperatures respectively 
(see §4.2l fc §4.31) . The eighth and tenth columns list the total column density for the assumed 
lower and upper rotation temperatures in cm -2 (see §4.31) . The fits to the spectra are chi 
squared gaussian fits and only those lines for which we have a full spectral profile have been 
fit. 

We did not set out to detect the second, narrowband, CH 3 OH transition but it was 
fortuitously at the edge of one of our spectral windows. Unfortunately we did not observe 
the entire line of this higher energy (J=19 5i i 5 — 20 4] i 6 ) CH 3 OH transition. We present a map 
of its emission at a wlsr of 5.5 km s -1 in Figure [3j The emission is restricted to the hot core 
region and shows a slight extension to the southwest, similar to that of the other CH3OH 
transition. 
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3.3. Dimethyl Ether [(CH 3 ) 2 0] 



(CH-Q oO is a typical l arge O-bearing species and has been associated with the compact 



ridge (e.g. iLiu et al 



20021 ). Figure UK shows the (CH 3 ) 2 map at a wlsr of 7.6 km s 



There clearly is some extended emission as indicated by the large negative contours. There 
is significant emission from IRc6 with a notable hole near Orion-BN, and in regions near 
IRc5. As with CH3OH there is no notable emission from the compact ridge its elf. Also there 



i s som e emission from the hot core-SW. This emission is similar to that seen by lBeuther et al. 



( 120051 ) at 865 /im, however they do not see the "C" shape around BN, but rather see more 
point like emission from IRc6. This may be explained by the factor of ~2 difference in E u 
between the observed transitions. Figures Hb and c show the (CH 3 ) 2 spectra from IRc6 
and IRc5, respectively, and Table [2] shows the fits to the spectra. The IRc6 emission is best 
fit by two gaussians, one narrow the other wide. The wide absorption feature, near km 
s _1 , seen in both spectra are due to the sidelobes of the C2H5CN (J=27o, 2 7 — 26o,26) emission 
from the hot core. 



3.4. Methyl Formate [HCOOCH3] 

Similar to (CH^O, HCOOCH3 has traditionally been associated with the compact 
ridge. Both the A and E torsional states of the HCOOCH3 J=20i i ig — 191,18 transition 
were detected. Figure [5^ shows the v lsr=7.6 km s -1 map of the A torsional state. Unlike 
(CH 3 ) 2 0, the strongest emission from HCOOCH3 comes from the IRc5 region and not IRc6. 
Additionally there is some emission from the hot core-SW and no detectable emission from 
the compact ridge. A s with (CHa)oO, ther e are some notable differences between our map 



and that observed by iBeuther et al.l (120051 ). especially the low level of emission we observe 
toward IRc6. This difference may be explained by the large (nearly factor of 3) difference 
in E u of the observed transitions or by our single channel maps vs their velocity-integrated 
maps. The peak in the northeast of our map is most likely from sidelobes that were not 
completely removed. Figures [5b and c show the HCOOCH3 spectra toward IRc6 and IRc5, 
respectively. The velocity scale is in reference to the A torsional state transition. The IRc5 
emission appears to have two peaks at ~7.6 and 9.3 km s _1 . 



3.5. Ethyl Cyanide [C 2 H 5 CN] 



C 2 H 5 CN is a t ypical large N-bea ri ng species and has been associated with the hot core 
for many years (e.g. iBlake et al.lll987l ). iBlake et al.l (119961 ) and lBeuther et al.l (120051 ) showed 
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an extension in C2H5CN emission toward IRc7. The map of C2H5CN 3=25^^2 — ^3,21 
emission, shown in Figure [6^, clearly resolves this second component for the first time. It 
also shows an extended ridge of emission running northeast-southwest along the hot core and 
in the direction of known outflows. Additionally, there is a peak of emission to the north of 
the hot core, a region which has no detected continuum emission. The emission from IRc7 
may be coming from two (or more) partially resolved sources. Figures [6b and c show the 
C 2 H 5 CN spectra toward IRc7 and the hot core, respectively. The C 2 H 5 CN spectra from IRc7 
were best fit by two gaussians, one narrow and one wide, while the hot core spectra were 
best fit by a single gaussian. The second line (t>LSR~-9.5 km s _1 ) in the hot core spectra 
may be a U-line (unidentified) or it may be a velocity component of C2H5CN. This line does 
not have a large spatial distribution and is confined to the southwest corner of the hot core. 
The other transition of C2H5CN we observed (J=27o,27 — 26o,26) maps nearly identically to 
this transition. 

One notable feature of the C2H5CN emission is a strong peak at a ^lsr of 12 km s -1 just 
south of IRc6. Figure [7^ shows the channel map of this emission and Figures [7b and c show 
the spectra from near the hot core and IRc6 respectively. There is no notable continuum 
associated with this emission and nothing else like it appears in any of the other transitions 
we observed. 

While there is no notable velocity gradient along the ridge of the hot core there are some 
velocity gradients across the hot core near its southwestern tip. Figure [H] shows a position- 
velocity (p-v) diagram of C2H5CN. The inset is a map of C2H5CN emission at a vlsr of 5 km 
s _1 (a zoom in of Figure [6^). The line denotes the p-v cut and the "+" denotes the center 
position of the cut (a(J2000) = 05 /l 35 m 14 s .4, <5(J2000) = -5 22'31".5). The abscissa of the 
p-v diagram is offset arcseconds from the center point and the ordinate is vlsr m km s _1 . 
The positions of the hot core and IRc7 are noted in the diagram as is the position of the U- 
line. The contours are the same as in Figure EJ The hot core contains not only C2H5CN gas 
at 5.0 km s _1 but some higher velocity gas, near a vlsr of 7.5 km s -1 (the typical compact 
ridge velocity). The diagram also shows that the hot core and IRc7 are connected by a ridge 
of lower velocity gas. This ridge also connects with the U-line indicating that it might be 
C2H5CN at a lowe r velocity. The U-li ne may also be the 1 z/ 6 J=25 — 24 le transition of HC3N 



at 227.7916 GHz flMiiller et all 120051 ). a rather high energy line (£ u =859 K), which could 



explain its compactness. 
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3.6. Acetone [(CH 3 ) 2 CO] 

(CH3)2CO, while a large O-bearing species, does not map like any of the other molecules 
we observed. Figure [9^, shows a map of the (CHa^CO emission. Two distinct emission 
features are seen. One at hot core-SW, the other near IRc7. Figures [9b and c show the 
spectra toward IRc7 and the hot core, respectively, and Table [2] lists the parameters of the 
spectral fits. All three spectral features (AE/EA, EE, & AA) are seen at both positions. 
Both spectra are best fit by two velocity components, ~5 km s _1 and ~8 km s -1 , typical 
vlsr values for the hot core and compact ridge. This is the largest known molecule that 
shows both velocity components, which are co-spatial at our resolution. Unlike C2H5CN 
there are no velocity gradients in any of the (CH 3 ) 2 CO emission. 

3.7. Formic Acid [HCOOH], Acetaldehyde [CH 3 CHO], and Ethanol [C 2 H 5 OH] 

The HCOOH, CH3CHO, and C2H5OH transitions listed in Tabled] were not detected 
above their respective 3 a detection limits. This may be due to these molecules having a 
generally extended emission and thus being resolved out by the array, or the transitions we 
observed were not favorable, or both. 

3.8. Lines Detected in Wideband Windows 

Several molecular transitions were strong enough to detect in our wideband windows 
(~500 MHz bandwidth, ~ 41.3 km s _1 resolution). Ten of the twelve transitions from 
six molecules detected in the wideband channels are listed in the lower part of Table [ifE 
Figure [10] shows the channel maps of the transitions. The 1 a rms noise level for all maps 
in this figure is 29 mJy bm -1 . Note that all maps in this figure are one channel wide, and 
cover all the velocity components of the transition. 

Figure [TUh shows the OCS J=19-18 emission. There are notable peaks and valleys all 
throughout the region indicating widespread features that are resolved out by the array, in 
addition to the smaller scale structures that are seen. For such a low temperature transition 
from a simple molecule this is not a surprising result. Figure [TUb shows the CH3OH J=102,g — 
93,6 emission. The emission is very similar to that of the CH 3 OH J=8_i g — 7qj transition. 

Figure [TUb shows the HCOOCH3 J=202,ig — 192, is A and E emission, which maps very 



2 The other two transitions are from HCOOCH3 and were also detected in the narrow spectral windows 
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similarly to the J=20i ) i9 — 19i,i8 A and E emission, with the exception that these HCOOCH3 
transitions have a strong peak near IRc6, which may indicate contamination from another 
species in the wideb and channel or may be due to velocity-averaging (similar to that of 



Beuther et all (|2005[ )). Figure OTOti-f show the l v n J=25 , 25 - 24 , 24 and J=24 2i23 - 23 2i22 
and 1 z/i5 J=25o, 2 5 — 24 , 2 4 transitions of vinyl cyanide [C 2 H 3 CN], respectively. These are the 
only transitions of C 2 HsCN detected and all are from vibrationally excited states. They map 
very similarly to the C 2 H 5 CN transitions, having strong emission from the hot core region 
and IRc7. One notable difference is that the C 2 H 3 CN emission shows two emission peaks 
in IRc7 while C 2 HsCN only hints at two spatially blended components. These two compo- 
nents have peaks at a(J2000) = 05 A 35 m 14 s .333, 5(J2000) = -05°22'31".5 and a(J2000) = 
05 fe 35 m 14 s .239, 5(J2000) = -05 o 22'31".l. 

Figure [TUg shows the J=26i )2 5 — 25i j2 4 transition of C 2 H 5 CN. While this transition maps 
similarly to the other detected C 2 HsCN transitions there is one significant difference. The 
emission near IRc6 seen here is not detected in any other C 2 HsCN transitions. This may be 
due to contamination from weaker lines of O-bearing species (several CH3CHO transitions 
fell within the same channel). Lastly, Figure [TOh shows the J=5-4 transitions of 13 CS. The 
map shows strong emission from the hot core, IRc7, and IRc6 and weaker emission from 
IRc5 and the hot core-SW. 



4. DISCUSSION 
4.1. Sources 

4.1.1. Hot Core and IRc7 

As can been seen from the maps of C 2 HsCN and C 2 HsCN, when observed under high 
spatial resolution, molecules that traditionally have been associated with the hot core have 
multiple strong components: one or two sources associated with IRc7 and a string of com- 
ponents associated with the hot core. 

The spectral line fits to the hot core features (listed in Table [2]) show a common ~5.2 
km s _1 velocity component (the canonical rest velocity for this source) for all species de- 
tected there. However, the two O-bearing species detected toward the hot core/hot core-SW, 
CH3OH and (CH 3 ) 2 CO, have a second, narrower, component near 7.5 km s _1 (a typical rest 
velocity for compact ridge species). How this second component is associated with the hot 
core is unclear, although since the CH3OH line is in absorption it is likely that this component 
is from gas along the line of sight. 
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The spectral line fits to the IRc7 features have fewer commonalities. The C2H5CN 
spectrum is best fit by a pair of gaussians at ~4 and ~5 km s _1 , both with distinctly 
different line widths. The O-bearing species in IRc7 also have 2 components, one near 5.1 
km s _1 the other near 7.5 km s _1 . The latter velocity component is in absorption in CH 3 OH, 
indicating that it is cooler and may be between us and IRc7. The C2H3CN maps give a strong 
indication of multiple cores which may have differing velocities. However they are unresolved 
in all other lines, so future observations will be needed to study these further. 



4-1-2. Compact Ridge, IRc6, and IRc5 

The spectral line fits to all the IRc6 features have a common component near 7.4 km 
s _1 . Additionally (CHa^O and CH3OH have secondary velocity components at ~9-10 km 
s _1 . As none of the CH3OH components are in absorption it is likely that all velocity 
components are warmer. Also it is evident in the maps that there are multiple unresolved 
components in IRc6, eac h of which may have a slightly different vlsr- Comparing our maps 



with those presented by iBeuther et al.l (120051 ) indicates that the structure that is observed 



in IRc6 may be more dependent on the E u of the transition observed, than it is for any of 
the other sources. This indicates that IRc6 may be at a notably different temperature, or 
have distinctly different physical conditions, than the other regions. 

As with IRc6, IRc5 has two distinct velocity components, one near 7.6 km s _1 , the other 
near 8.5 km s _1 , for all transitions observed toward this source. The CH3OH transition is 
in absorption, indicating that it may be between IRc5 and us. As with our other sources, 
IRc5 appears to have at least 2 partially resolved components, which may be the sources 
of the different velocity components. As no transitions of larger molecules were detected 
toward the compact ridge it is likely that many molecules typically associated with the 
compact ridge, via their i>lsr, ma y actually be in IRc5 instead. One exception that we know 
of is HCOOH. Obs ervations comparing the spatial distribution of (CH 3 ) 2 and HCOOH 



by iLiu et al.l (120021 ) show that HCOOH may be associated with the compact ridge as its 
emission is to the southeast of that of (CH 3 ) 2 and appears to be on top of the compact 
ridge. 



4.2. Rotation Temperature Estimation 



As we did not detect enough transitions of any one molecule to construct a reason- 
able rotation-temperature diagram we will assume a range of temperatures for the different 
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molecular species. From the 3mm spectral line survey of Friedel et al., in preparation, we 
will use temperatures of 100-150 K for O-bearing (compact ridge) species and 200-500 K for 
N-bearing (hot core) species. Note that for (CH 3 ) 2 CO we assume a temperature range of 
the other O-bearing species as its original detection gave rotation temperature of ~176 K, 
notably below the range from the hot core. 



4.3. Opacity and Column Densities 



High optical depth ( r > 1) can greatly reduce the observed intensity of a spectral line. 
From iRohlfs and Wilsonl (120001 ) the opacity at the line center can be calculated from 



-In 



T, 



MB 



J(T r ) - J (2.73 K) - T< 



c 



(1) 



where T MB = T^/r] M is the main beam brightness temperature of the line and rf M is the 
corrected main beam efficiency J(T) is defined as (IRohlfs and Wilsonl 120001 ) 



J(T) 



hv 



k e hv l kT - 1 



K. 



(2) 



J(2.73) is the cosmic microwave background, and Tq is the continuum brightness tempera- 
ture. For observations not done in a temperature scale one can convert from Jy beam -1 to 



K using (IRohlfs and Wilsonl l2000h 



T, 



MB 



1.22.r , 

% x 10 6 K, 

B0 a b v* 



(3) 



where Iq the is peak intensity of the line in Jy beam -1 , B is the beam filling factor (assumes 
to be 0.5, filling the beam, for all calculations here), 6 a and 6 b are the FWHM Gaussian 
synthesized beam axes in arcseconds, and v is the transition frequency in GHz. For emission 
lines the optical depth correction factor, C T can then be calculated ([Goldsmith fc Langer 
1999h : 

°r = 7^7- (4) 
1 — e T 

This factor is then applied to the column density equation (jSJ). The opacities for the detected 
transitions are given in the seventh (low temperature limit, 150 K for hot core/IRc7 and 100 
K for all others) and ninth columns (high temperature limit, 500 K for hot core/IRc7 and 
150 K for all others) of Table H 

With the assumption of LTE, to calculate the tot al beam averaged c olumn density 



((Nt)) from a single transition for an array, we have (e.g.. iFriedel et al.ll2005l ; iRemijan et al. 
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2005|) 



and 



Nu 
9u 



2MWC T 



T 



B6 a 6 b S^ T r - T b 



x 10 



20 cm" 2 



(N T ) 



^% TV c 

9u 



(5) 



(6) 



Here W is the integ rated line intensity i n Jy beam -1 km s -1 , Q r „ is t he rotational- vibrational 
partition function ( IFriedel et al.l 120051 : IWidicus Weaver et al.ll2005l ). E u is the upper state 
energy of the transition in K, T r is the rotational temperature in K, Sfi 2 is the product 
of the line strength and the square of the relevant dipole moment in Debye 2 , N u /g u is the 
upper state column density divided by its degeneracy (2J+1), and TJ, is the total background 
continuum temperature. The column densities for each transition are given in the eighth 
(low temperature limit, 150 K for hot core/IRc7 and 100 K for all others) and tenth (high 
temperature limit, 500 K for hot core/IRc7 and 150 K for all others) columns of Table [2j 

From Table [2] it can be see that most of our lines are fairly optically thin at our low 
temperature limit. However, several transitions, notably those of O-bearing species, are 
optically thick. At our high temperature limit it can be seen that more transitions are 
optically thin when compared to the low temperature limit. In general the column densities 
for all transitions are notably higher than those reported previously in the literature. The 
reasons for this are 1) we are resolving the sources and thus do not suffer from beam dilution 
and 2) we are accounting for the vibrational part of the partition function which can increase 
(Nt) by factors of 5 or more. 



4.4. Implications for Astrochemistry 



One of the big questions about the Orion-KL region is what is causing the O-N dif- 
ferentiation. One could argue that the O-N differentiation seen in Orion-KL is due to the 
differing dipole moments (and hence differing critical densities) of O and N bearing molecules 
(in general N bearing molecules have notably larger dipole moments). In order to investigate 
this possibility we calculate the critical density for the C2H5CN a nd (CH.Q^O tr ansitions we 



detected. The critical densities, be calculated as follows (ITielensI 



2005J): 



n. 



ul 



(7) 



where A u i is the Einstein A of the transitions and 7^ is the collisional de-excitation rate 
coefficient. Since the quantum levels of the C2H5CN and (CH 3 ) 2 transitions are dominated 
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by single radiative routes in and out, we can consider this to be a two level system and the 
summation collapses to single expressions. A u[ is defined as: 

2,hc 6 g u 

and we find for the observed transitions of C2H5CN (J=27o, 2 7— 26o,26) and (CH 3 ) 2 (J=13o,i 3 — 
12t 12 ), Ai./=9-8x l0~ 4 and 9.1xl0~ 5 s -1 respectively. The de-excitation coefficient is defined 



as (ITielensI 120051 ): 



'"' ^ ( 2^) 3/2 I °Uvye-^dv, (9) 



where /1 is the reduced mass of the system and a u i(v) is the collisional de-excitation cross 
section at the relative velocity v. If we assume the geometrical cross section^ of each molecule 
(a ul (v) « 1.3 xlO^ 15 cm 2 for C 2 H 5 CN andcx tti (t;) « 7.1 xlO^ 16 cm 2 for (CH 3 ) 2 0), the primary 
collision partner to be H 2 , and temperatures of 200 K and 100 K0for C 2 H 5 CN and (CH 3 ) 2 
respectively, then we find 7 U ; = 1.9 x 10~ 10 cm 3 s _1 for C 2 HsCN and jui = 7.4 x 10~ n cm 3 s _1 
for (CH 3 ) 2 0. Using equation ([7]) gives critical densities of 5.3 xlO 6 cm -3 for C 2 H 5 CN and 
1.2xl0 6 cm -3 for (CH 3 ) 2 0. This slight difference in critical densities is unlikely to be the 
cause of the differentiation. 

Others have argued that the differentiation could be due to differences in the age of the 
regions, the O-bearing species are liberated off grain surfaces by shocks or form in gas phase 
reactions in s hocked regions, w hile the N-bearing species take longer to form in gas-grain 



reactions (e.g. iBlake et al.lll987l ). This argument fits well with what is suspected about the 
region: outflows from the older (N-bearing) hot core are impacting the ambient gas and 
dust forming shocked regions where larger O-bearing species are most populous. However, 
our observations of (CH 3 ) 2 CO show that this O-N differentiation is not as simple as it may 
seem. Figure [TT] shows C 2 H 5 CN in grey (red in electronic edition) contours, (CH 3 ) 2 as 
dashed (green in electronic edition) contours, and (CH 3 ) 2 CO as black contours, overlayed 
(note that negative contours have been removed to reduce confusion). It shows that the 
(CH 3 ) 2 CO emission comes only from regions where both (CH 3 ) 2 and C 2 H5CN overlap. 
If we assume that the age difference argument is correct then it would seem that it takes 
significant processing (either on grain surfaces or in the gas phase) to form (CH 3 ) 2 CO, but 



3 This assumption is valid since the primary collision partner, H2, has no permanent dipole moment and all 
other interactions are very short range, thus the collisional cross section will be dominated by the geometrical 
cross section. 

4 Changing the assumed temperature will not affect the outcome much as n cr will go as T -1 / 2 (higher 
temperature gives a lower n cr ). 



-13- 



that it is also short lived as it seems to not be present in older regions where N-bearing 
molecules dominate. Recent models by Garrod, Widicus Weaver, & Herbst, in preparation, 
indicate that this may be true: (CH 3 ) 2 CO forms where the warm-up time scale is longer, 
indicating grain surface formation from reactions with secondary products. 

5. SUMMARY 

We have conducted high resolution A=lmm observations of the Orion-KL region with 
the CARMA array. Our observations studied the continuum of the region as well as eleven 
molecular species. Of these eleven species, eight were detected and we have presented their 
maps and spectra. These maps show that this region more complex than previously thought. 
The large N-bearing species (C2H5CN and C2H3CN) originate from multiple cores near the 
Orion hot core and IRc7, while large O-bearing species ((CH 3 ) 2 and HCOOCH3) originate 
near IRc5 and IRc6, and not the compact ridge as had been assumed before. Additionally 
we see that (CH 3 ) 2 CO, unlike other large O-bearing species, does not originate from IRc5 
or IRc6, but rather from the hot core-SW and IRc7, the only regions where both large 
N-bearing and O-bearing species are present. 

We would like to thank S. Widicus Weaver and L. Looney for many helpful comments 
and discussions. We also thank an anonymous referee for a very favorable review of this 
work. This work was partially funded by NSF grant AST-0540459 and the University of 
Illinois. Support for CARMA construction was derived from the states of Illinois, California, 
and Maryland, the Gordon and Betty Moore Foundation, the Kenneth T. and Eileen L. 
Norris Foundation, the Associates of the California Institute of Technology, and the National 
Science Foundation. Ongoing CARMA development and operations are supported by the 
National Science Foundation under a cooperative agreement, and by the CARMA partner 
universities. 
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Table 1. Molecular Parameters of Observed Lines 



Quantum Frequency E u S/i 2 a X &b RMS Noise Missing 

Molecule Numbers (MHz) (K) (D 2 ) (" X ") (mjy bm" 1 ) Flux (%) a Ref. 



Primary targets in narrowband windows 



H000CH3 


20i,i 9 


- 19i,is A 


226,778.764 


(10) 


120 


51.4 


2.49 x 0.90 


185 


~ 15 


1 




20i,ig 


- 19i,i8 E 


226,773.152 


(9) 


120 


51.4 


2.49 x 0.90 


185 


~ 15 


1 


(CH 3 ) 2 


13o,13 


- 12i i:L2 AA b 


231,987.779 


(3) 


81 


16.9 


2.42 x 0.82 


252 


~ 10 


2 




13o,13 


— 121,12 EE b 


231,987.856 


(3) 


81 


16.9 


2.42 x 0.82 


252 


~ 10 


2 




13o,13 


— 121,12 AE/EA 


231,987.933 


(3) 


81 


16.9 


2.42 X 0.82 


252 


~ 10 


2 


C2H5CN 


27o,27 


— 26o,26 


231,990.409 


(3) 


157 


399.4 


2.50 x 0.83 


252 


< 5 C 


3 




25 3j2 2 


- 24 3 ,2i 


227,780.971 


(2) 


150 


365.4 


2.42 x 0.83 


1 79 


< 5 C 


?> 


(CH 3 )2CO 


23o,23 


- 22i,22 AE Q 


230,170.265 


(20) 


135 


22.4 


2.47 x 0.82 


196 




1 




23(1,23 


- 22i,22 EA d 


230,170.313 


(17) 


135 


22.4 


2.47 x 0.82 


196 


c 


4 




23o,23 


- 22i,22 EE 


230,176.727 


(14) 


135 


22.4 


2.47 x 0.82 


196 


c 


4 




23o,23 


- 22i,22 AA 


230,183.083 


(21) 


135 


22.4 


2.47 x 0.82 


196 


c 


4 


CH3OH 


8-1,8 ' 


- 7o,7 


229,758.811 


(15) 


89 


5.0 


2.45 x 0.88 


268 


~ 90 


5 




19 5 ,15 


— 204,16 


229,864.221 


(46) 


578 


5.7 


2.43 x 0.88 


188 


f 


5 


HCOOH 


10l,9 - 


-9i, 8 


231,505.705 


(23) 


64 


19.3 


2.38 x 0.89 


279 


g 


3 


CH3CHO 


111, 11 


- 10o,io 


229,775.029 


(11) 


61 


8.4 


2.45 x 0.88 


284 


g 


6 


C2H5OH 


10 2 , 9 - 


-9i,s 


226,661.701 


(100) 


51 


17.0 


2.47 x 0.90 


218 


g 


3 








Transitions detected 


in wideband windows 








ocs 


19-18 




231,060.983 


(10) 


110 


9.7 


2.44 x 0.83 


29 


h 


3 


CH30H 


10 2 ,9 - 


-93,6 


231,281.150 


(15) 


165 


2.67 


2.44 x 0.83 


29 


h 


5 


HCOOCH3 


20 2 ,i9 


- 192,18 A 


226,713.028 


(9) 


120 


51.4 


2.48 x 0.84 


29 


h 


1 




20 2 ,i9 


- 192,18 E 


226,718.748 


(10) 


120 


51.4 


2.48 x 0.84 


29 


h 


1 


C2H3CN Vi 


242,23 


— 232,22 


226,835.668 


(16) 


488 


346.7 


2.48 x 0.84 


29 


h 


3 




24i,23 


— 23l,22 


231,101.164 


(14) 


486 


348.3 


2.44 x 0.83 


29 


h 


3 




25o,25 


- 24 ,24 


231,145.707 


(14) 


190 


363.1 


2.44 x 0.83 


29 


h 


3 


Vl5 


25o,25 


— 24o,24 


231,235.325 


(75) 


635 


363.3 


2.44 x 0.83 


29 


h 


3 


13 cs 


5-4 




231,220.996 


(200) 


33 


19.3 


2.44 x 0.83 


29 


h 


3 


C 2 H 5 CN 


26i,25 


-25i,24 


231,310.420 


(50) 


153 


383.0 


2.44 x 0.83 


29 


h 


3 



a When compared with the results from lsutton et all (|l985h . 
b These (CH 3 ) 2 lines are blended. 

c Calculations indicate that little or no flux is missing, but since the features are unresolved in our observations and source sizes 
cannot be determined the values are upper limits. 

d These (CH 3 ) 2 CO lines are blended. 

°Not detected bv lSutton et all d 19851) . 

f The full line was not detected in our observations, so no comparison could be made. 
g This line was not detected in our observations, so no comparison could be made. 
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h Since this line was spectrally unresolved, no comparison can be made. 



References. — (1) lOesterling et al.l l|l99Sh : (2) IGroner et all il998ll ; (3) IPickett et all <1998h : (4) IGroner et all j2002h ; (5) 
IXu fc Lovasl Jl99Ttl : f6) lLovas I 120041 1 
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Table 2. Fits and Column Densities of Observed Lines 



Low T b High T c 











^LSR a 






Aim 






IVIoIccliIc 


Source 


f lg. 


(jy Dm j 


(km s ) 


(km s ) 


T 


(cm -2 ) 


T 


(cm~ 2 ) 


C 2 H 5 CN 


IRc7 




2.5(2) 


4.9(1) 


2.8(2) 


0.36(4) 


3.2(4)xl0 16 


0.13(1) 


1.9(3) xlO 17 




IRc7 


■Bp 


0.8(2) 


3.9(8) 


11.6(11) 


0.10(3) 


3.7(10)xl0 16 


0.04(1) 


2.4(6) xlO 17 




Hot Core 




3.9(1) 


5.1(2) 


6.3(1) 


0.74(4) 


1.5(l)xl0 17 


0.21(1) 


7.0(3)xl0 17 


(CH 3 ) 2 


IRc6 


Wp 


1.8(4) d 


7.1(1) 


1.5(4) 


0.63(20) 


9.0(3) xlO 17 


0.35(10) 


1.3(5)xl0 18 




IRc6 




2.8(3) 


9.1(4) 


6.0(6) 


1.29(28) 


7.4(13)xl0 18 


0.62(9) 


9.3(14)xl0 18 




IRc5 


it 


3.4(4) 


7.6(1) 


1.1(2) 


1.77(65) 


1.9(6)xl0 18 


0.78(15) 


2.1(5)xl0 18 




IRc5 




0.9(3) 


8.1(6) 


2.5(14) 


0.25(10) 


6.0(40) xlO 17 


0.16(6) 


9.6(64) xlO 17 


HCOOCH 3 e 


IRc6 




0.5(1) 


7.7(4) 


3.3(8) 


0.13(3) 


3.1(10)xl0 16 


0.08(2) 


4.6(15)xl0 16 




IRc5 


Et 


3.7(3) 


7.5(0) 


1.3(1) 


2.11(7) 


2.0(5) xlO 17 


0.83(11) 


1.8(2)xl0 17 




IRc5 


Et 


2.4(4) 


9.3(0) 


0.9(1) 


0.84(23) 


3.6(8) xlO 17 


0.45(10) 


4.4(9) xlO 17 


CH3OH 


IRc5 




4.1(5) 


7.8(2) 


6.6(6) 


4.37(1030) 


1.5(33) xlO 18 


1.00(22) 


9.2(17)xl0 17 




IRc5 


lit 


-4.1(3) 


8.6(1) 


1.4(3) 


f 




f 






IRc6 




2.3(4) 


7.1(2) 


1.4(7) 


0.80(22) 


5.9(32) xlO 16 


0.44(10) 


8.6(46) xlO 16 




IRc6 




5.0(4) 


8.7(2) 


0.6(3) 


g 




1.47(29) 


1.2(6)xl0 17 




IRc6 




1.5(3) 


10.3(2) 


1.3(4) 


0.45(12) 


3.1(ll)xl0 16 


0.26(6) 


4.8(18)xl0 16 




Hot Core 


2.1 


2.8(2) 


5.6(2) 


5.5(4) 


0.42(4) 


7.4(8) xlO 17 


0.13(1) 


4.7(5) xlO 18 




Hot Core 


2.1 


-1.8(4) 


7.4(1) 


1.1(3) 


f 




f 






IRc7 


m 


1.3(3) 


4.9(8) 


6.9(16) 


0.16(4) 


3.5(ll)xl0 17 


0.06(1) 


2.5(9) xlO 18 




IRc7 


m 


-1.2(6) 


7.0(3) 


1.2(9) 


f 




f 




(CH 3 ) 2 CO 


IRc7 


9 3 


1.3(2) h 


8.0(2) 


2.4(4) 


0.37(6) 


3.2(7)xl0 17 


0.23(3) 


5.5(12)xl0 17 




IRc7 




0.7(2) h 


5.3(4) 


2.4(10) 


0.18(5) 


1.6(8)xl0 17 


0.12(3) 


2.8(14) xlO 17 




Hot Core-SW 





1.7(2) h 


7.7(2) 


2.1(4) 


0.75(15) 


5.7(13)xl0 17 


0.38(6) 


8.1(18)xl0 17 




Hot Core-SW 




0.4(2) h 


5.3(8) 


2.5(6) 


0.14(7) 


1.3(7)xl0 17 


0.08(4) 


2.1(ll)xl0 17 


a All uncertainties listed are 


2(7, an 


uncertainty of 


means the there are not enough sig 


nificant figures to display its true value. 



b ow temperature for hot core/IRc7 is 150K, for all other sources is 100 K 



c High temperature for hot core/IRc7 is 500K, for all other sources is 150 K 
d /o for blend of all torsional states. 

°Quantities are for both A and E torsional states, as they are assumed to have equal populations. 
'Opacity and column density not calculated for absorption components. 

g Opacity calculation gave negative values, likely due to too low an assumed rotation temperature. 

h 7o for EE torsional state only; for other states as follows: 7 (AA) = 0.5io(EE), / (EA) = 7 (AE) = 0.25/ (EE). 
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Fig. 1. — Continuum of Orion-KL at 231 GHz. Known infrared and radio sources are denoted 
with "+" and labeled. Contour levels are ±3er, ±5<r, etc., where the 1 a rms noise level is 
11.3 mJy bm _1 . The synthesized beam is in the lower left corner. 
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Fig. 2. — Map and spectra of CH 3 OH. a) Map of the J=8_i ;8 — 7 j, fLSR=5-5 km s -1 emission 
overlayed on gray-scale continuum. Contours are ±3cr, ±5<r, ... (a = 0.268 Jy bm _1 ). The 
synthesized beam is in the lower left corner of the map and "+" marks denote the same 
sources as in the continuum map (see Figured]), b) Spectra toward IRc5. Dashed lines mark 
vlsr=7.6 and 9.6 km s -1 . "I" bar denotes 1 a rms noise level, c) Spectra toward IRc6. d) 
Spectra toward hot core, e) Spectra toward IRc7. 



20 -10 10 20 



-21 - 



C W 19 5,19- 20 4,16 




+ 



14 s 5 14 s 
a (2000) 



13.5 



Fig. 3. — Map of CH 3 OH J=19 5j i5 — 20 4i i 6 , w lsr=5.4 km s -1 emission overlayed on gray-scale 
continuum. Contours are ±3cr, ±5cr, ... (a = 0.188 Jy bm _1 ). The synthesized beam is in 
the lower left corner of the map and "+" marks denote the same sources as in the continuum 
map (see Figured]). 



-22 - 



20 



25 



30 



o 35 



40 - 



45 



-05 22 50 



a) (CH 3 ) 2 13 0i13 -12 
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Fig. 4. — Map and spectra of (CH 3 ) 2 0. a) Map of the t>LSR=7.6 km s -1 emission overlayed 
on gray-scale continuum. Contours are ±3cr, ±5<r, ... (a = 0.252 Jy bin -1 ). The synthesized 
beam is in the lower left corner of the map and "+" marks denote the same sources as in the 
continuum map (see Figured]), b) Spectra toward IRc6. Dashed lines mark 5.0 and 7.6 km 
s _1 (typical rest velocities for this region). "I" bar denotes 1 a rms noise level, c) Spectra 
toward compact ridge. 
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a (2000) 



13.5 




Fig. 5. — Map and spectra of HC00CH 3 . a) Map of the A torsional state, t>LSR = 7.6 km 
s _1 emission overlayed on gray-scale continuum. Contours are ±3<r, ±5<r, ... (a = 0.185 
Jy bm _1 ). The synthesized beam is in the lower left corner of the map and "+" marks 
denote the same sources as in the continuum map (see Figure [T]). The emission component 
to the northeast of the hot core may be a sidelobe that was not completely removed, b) 
Spectra toward IRc6. Dashed lines mark i>lsr = 7-6 and 9.2 km s _1 . "I" bar denotes 1 a rms 
noise level. The second set of spectral features near 15 km s -1 are the emission from the E 
torsional state of this transition, c) Spectra toward IRc5. 
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Fig. 6. — Map and spectra of C 2 H 5 CN. a) Map of the J=25 3>2 2 — 24 3j2 i, fLSR=5.0 km s _1 
emission overlayed on gray-scale continuum. Contours are ±3<r, ±5<r, ... (a — 0.179 Jy 
bm _1 ). The synthesized beam is in the lower left corner of the map and "+" marks denote 
the same sources as in the continuum map (see Figured]), b) Spectra toward IRc7. Dashed 
lines mark vlsr = 5.0 and 7.6 km s" 1 (typical rest velocities for this region). "I" bar denotes 
1 a rms noise level, c) Spectra toward the hot core. 
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Fig. 7. — Map of C2H5CN at a wlsr of 12 km s _1 . a) Contours are the same as in Figure El 
Note the emission south of IRc6. b) Spectra toward the emission peak south of the hot core, 
c) Spectra toward the emission peak near IRc6. 
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-5 5 

Offset (") 

Fig. 8. — Position- velocity diagram of C2H5CN. The inset is a zoom in on the hot core/IRc7 
emission, at a i>lsr of 5 km s _1 , shown in Figure [61 The line denotes the p-v cut and the "+" 
denotes its center (a(J2000) = 05 fe 35 m 14 s .4, <5(J2000) = -5°22'31".5). The abscissa of the 
diagram is offset arcseconds from the p-v center and the ordinate is ^lsr m km s _1 . Labels 
indicate the position of the hot core and IRc7 along with the potential U-line. The contours 
are the same as those in Figure [6] 
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Fig. 9.— Map and spectra of (CH 3 ) 2 CO. a) Map of the J=23 , 2 3 - 22i j22 EE, w L sr=5.0 km 
s _1 emission overlayed on gray-scale continuum. Contours are ±3(7, ±5cr, ... (a = 0.196 Jy 
bin -1 ). The synthesized beam is in the lower left corner of the map and "+" marks denote 
the same sources as in the continuum map (see Figured]), b) Spectra toward IRc7. Dashed 
lines mark i>lsr = 5-0 and 7.6 km s _1 (typical rest velocities for this region). "I" bar denotes 
1 cr rms noise level c) Spectra toward hot core. 
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a (2000) 

Fig. 10. — Maps from the transitions detected in the wideband windows overlayed on gray- 
scale continuum. Contours are ±3cr, ±5(7, ... (1 o = 29 mJy bin -1 ). The synthesized beam 
is in the lower left corner of the map and "+" marks denote the same sources as in the 
continuum map (see Figured]). 



-29- 




05 h 35 m 15 s 14 s 5 1 4 S 1 3 S 5 

a (2000) 



Fig. 11. — Overlayed maps of C2H5CN in grey (red in electronic edition) contours, (CH3) 2 
in dashed (green in electronic edition) contours, and (CH 3 ) 2 CO in black contours. Note that 
negative contours have been removed to reduce confusion. The (CH 3 ) 2 CO emission only 
appears where both C 2 H 5 CN and (CH 3 ) 2 emission overlap. 



